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Abstract: We use effective field theory tools to study non-conservative effects in the 
gravitational 2-body problem in general spacetime dimension. Using the classical version 
of the Closed Time Path formalism, we treat both the radiative gravitational field and 
its dynamical sources within a single action principle. New results include the radiation- 
reaction effective action in arbitrary dimensions to leading and +1PN orders, as well as 
the generalized Quadrupole formula to order +1PN. 
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1 Introduction 

In [1] an effective field theory (EFT) formalism for the simultaneous, action-level treatment 
of radiation and radiation-reaction (RR) effects was developed (see also [2] for a pedagogical 
introduction). The EFT was explicitly constructed for systems of localized objects coupled 
to scalar, electromagnetic (EM) and gravitational fields. The method was generalized 
to arbitrary spacetime dimensions in the scalar and EM cases [3]. In particular, for a 
single point particle, it yields the higher-dimensional analogue of the Abraham-Lorentz- 
Dirac (ALD) Self-Force (SF). In this article we complete the generalization for General 
Relativistic (GR) post-Newtonian (PN) systems in arbitrary dimension by applying our 
method and constructing the associated EFT. The GR case is more involved than the 
(free) scalar and EM cases in two main respects: first, the higher spin of the gravitational 
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field which complicates the construction of gauge-invariant fields and corresponding sources, 
both central features of our formalism; and second, the intrinsic non-linearity of GR. 

The PN approximation of the gravitational 2-body problem in 4 spacetime dimensions 
has been studied, over the past few decades, to very high accuracy (for a wide review of the 
field see [4]). The state of the art is the +4PN (corrected to order v 8 ) effective action, re¬ 
cently completed in [5]. Since Goldberger and Rothstein’s groundbreaking paper [6] much 
progress was made on the EFT approach to the PN binary problem (including [7-14]); for a 
more detailed review section see [1]. Less attention has been given to gravitational radiation 
and radiation-reaction in spacetime dimensions other than 4 1 ; these will be the focus of this 
paper. Although not directly relevant for astrophysical gravitating binaries, we feel that 
studying this rather fundamental problem is well-motivated both because an understanding 
of a system’s behaviour in general dimension complements and enhances its understanding 
in 4 d (see for example [19-21]), and because higher-dimensional gravitational scenarios 
frequently emerge in theoretical physics, for example in string theory. Gravitational radi¬ 
ation in higher dimensions was treated in [22-24], including in EFT methods. Specifically, 
the leading-order expression for gravitational-wave emission in higher dimensions, i.e. the 
equivalent of Einstein’s Quadrupole formula [25], was computed. Eq. (4.25) of [24] (also 
(38) & (39) of [22]) gives the energy output in gravitational quadrupole radiation in d 
dimensions as 


dE 

du 


2 2-«* 7r -(d-5)/2 d (d-3) d+2 
d (d-2)(d+l)T [*±] 


\Qij(u )\ 2 


( 1 . 1 ) 


In this paper we develop an effective action for the radiation reaction of gravitation in 
general dimension (3.2). From this action we derive both the self-force and the dissipated 
power, up to subleading (+1PN) order. At the leading order, the energy output we find 
(3.6, 3. 7) matches exactly the result (1.1). As an additional test of our results at +1PN, 
we substituted d = 4 and compared with [1, 26], finding a match. 

It is important to stress a relevant qualitative difference between d = 4 and d > 4 
spacetime dimensions. For PN astrophysical binaries RR is weak, hence it must influence 
the system for a long time for its effect to be significant; this is the case for bound orbits. In 
higher dimensions, it is well known that there are no stable bound gravitational orbits for 
a binary system of gravitating compact objects 2 . Nevertheless, high-dimensional systems 
can be stabilized by non-gravitational forces (in particular short-range repulsive forces); 
RR may also have a substantial effect for trajectories close to the unstable circular orbit, 
as well as for systems with more than 2 bodies. As we are neither specifying the sources’ 
trajectories nor trying to solve for them, the question of binding (and its mechanisms) do 
not affect the results of this paper. 


1.1 Method 

In this work we treat the gravitational 2-body problem in arbitrary spacetime dimension, 
and calculate explicitly its radiation source multipoles, the RR effective action and im¬ 
portant physical quantities derivable from it: the outgoing radiation, dissipated energy 

1 For earlier work on a fully relativistic treatment of SF for fields of various spin, see [15-18]. 

2 For point particles, for example, there exist unstable circular orbits. 
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Figure 1. The two relevant zones: the system zone is on the left, with a typical stationary-like 
field configuration. The radiation zone with its typical out-spiraling waves is on the right. 

and RR force acting on the system. Our method [1] divides the problem into 2 zones, 
the system zone and the radiation zone (fig. 1) , each with different enhanced symmetries. 
The system zone enjoys approximate stationarity (time independence) since by assumption 
all velocities are non-relativistic; we thus use the “Non-Relativistic Gravitational” (NRG) 
fields [11] to describe it. The radiation zone enjoys an approximate spherical symmetry, 
as from its point of view the system has shrunk to a point. Hence we use gauge-invariant 
spherical field variables [27-31] as in [1] and unlike the plane-wave decomposition used 
in most previous EFT works [7-10, 14]). An important ingredient of our method is the 
matching of these two (system and radiation) zones at the level of the action, using ‘2-way 
multipoles’ [1]. These are degrees of freedom we introduce - ‘integrate in’ - to couple the 
two zones. From the radiation zone point of view we think of them as sources situated at 
the origin while in the system zone view they reside at infinity. At the end of the day we 
integrate out (eliminate) all the other degrees of freedom in the problem and remain with 
an effective action which is a function of these multipoles precisely. Schwinger-Keldysh 
field doubling [32] , which was beautifully adapted to the classical context in general in [33] 
and was used in studies of the binary problem in [1, 14], is an essential ingredient since it 
allows the derivation of dissipative effects from action principles. As we have shown in [3] 
and is further exemplified here, a clear advantage of our formalism is its ability to naturally 
extend to general spacetime dimensions. 

1.2 Conventions and nomenclature 

We follow the conventions of [3]. Thus the flat d-dimensional spacetime metric r/^’s sig¬ 
nature is mostly plus, and also D := d — 1, d := d — 3. Regarding the d+ 1 dimensional 
unit sphere we define to be its area; on it goo' is the metric, Dq is the covariant 

derivative 3 , and Aj +1 = DqD n is the Laplace-Beltrami operator. The eigenfunctions of 

3 At some points we will be interested in the covariant derivative with respect to an angular variable £7 
in the full d-dimensional space, and we will mark it Dq. Dn will be reserved for the more common case, 
of variation restricted to the , 1 sphere. Likewise the full d-dimensional metric, when it relates spherical 
coordinates, will be marked g^/- We note that g = r 2 gctn', where r is the radial coordinate. 
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this operator on the d+1 unit sphere will be given by various multipoles enumerated by an 
order l\ in treating their eigenvalues we shall make use of c s = £(£ + d) and c s = c s — (d+1). 
We shall designate by e the sector (scalar, vector, tensor) of different multipoles, with either 
e G {S, V, T} or e E {0,1, 2} respectively. 

Lower case Greek letters stand for spacetime indices {0..D}, lower case Latin letters 
for spatial indices {1..D}, upper case Greek letters for indices on the sphere {l..(d+ 1)}, 
Hebrew letters (M,H) for different vectorial and tensorial multipoles on the sphere, and 
upper case Latin letters for spatial multi-indices, i.e. I = Ig := (%£... ig) where each 
ik £ {1..D} is an ordinary spatial index, and £ is the number of indices. We use the 
multi-index summation convention with implied 

Pi Qi ■= Y P h := Y Qh-ii > ( L2 ) 

i I ’ 

as well as the multi-index delta function defined through dj e j e := £\5i 1 j 1 .. ■Si e j i , so that 
factors of £! are accounted for automatically. 

While as usual c = 1, we specifically wish to keep the gravitational constant in d 
dimensions, marked Gd- We choose it and the normalization for the gravitational action 
so that the Newtonian potential is always —GdM/r d , and the action is given by (2.1). It 
is related to the commonly used constant G^ 4 by Gd = ^ — G^. Both definitions 

identify for d = 4 [24, 34]. The constant Gd also has units of L d /M. 

We use the convention that the Feynman rules are real [35], as befits a non-quantum 
field theory. We also denote 

(L3) 


2 From Einstein’s Action to Feynman Rules 


2.1 Action for metric perturbations 

In GR, the action for the metric g^ u coupled to matter is given by the sum of the Einstein- 
Hilbert (EH) Action and a matter term: 


S 



—^- R + C-m i 

2(d+l)Q^ +1 Gd 


( 2 . 1 ) 


where Cm is the matter Lagrangian. We shall calculate the generation and reaction of 
Gravitational Waves (GWs) propagating out to infinity in asymptotically flat space, in 
spherical coordinates. The metric splits to g fllJ = g fllJ + h^, where g^ describes the 
background flat spacetime in spherical coordinates 5 and is a perturbation on it. The 
EH term is the kinetic part of the action, and describes propagation of GWs. The matter 

4 In terms of which the overall prefactor of the Einstein-Hilbert action is 167r 1 Gd in arbitrary dimension 
and black-hole entropy is given by 

5 See appendix A; henceforth, covariant derivatives are associated with this metric. 
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term describes the generation of GWs from (matter) sources. We expand the kinetic 
(source) action to quadratic (linear) order in h^ v to find 


Seh 

Smat 


1 

4(d+l)n^ +1 Grf 



1 

2 


h a 1 h P 

n [or, p];7 


+ hfj? h 


a 

7 ];a ’ 



x^gh^ . 


( 2 . 2 ) 

(2.3) 


We use the action composed of (2.2) and (2.3) to construct our far zone action in terms 
of spherical variables. As we compute corrections only up to +1PN, and as far zone 
nonlinearities enter only at orders higher than +1PN (see discussion in section 3.3), we 
need not take them into account here; this justifies expanding the action to quadratic 
order in h. The near zone (a.k.a “system zone”) nonlinearities do enter at +1PN, and are 
treated in section 2.3. 

We use spherical harmonics to decompose the fields and sources to 11 families of fields, 
comprised of the 7 scalar families htt, h tr , h rr , hts, h r $, hs, hs , the 3 vector families h t $, 
h r x, h$, and the single tensor family /ijq, defined through (see appendix B) 


htt 



htr — 


h^n^e 


—iut 


hr-r - 


hl^n^e 


—iut 


htsi = ^2 {ht^danL + hftfnkn) e luJt , Kn {h^dan L + h%n% a ) e lut 


hnn< = 


hsfnnsy + a , + h^ u n^ nn , + h^n^nn' 


\Lu~L 


Lu^L 


^Lu^L 


—iut 


(2.4) 


where we use the scalar multipoles n L , the vector multipoles and the tensor multipoles 
n m , nn 1 described in appendix B. They are dimensionless and depend only on the angular 
coordinates. We shall at times omit L, oj indices for brevity. 

Substituting the new fields (2.4) into the homogenous action (2.2) and using the def¬ 
initions (B.4), derivative relations (B.6-B.11) and normalization relations (B.12) we find 
the homogenous action decomposes into 3 independent sectors, namely spherical scalar, 
vector and tensor. Following [27] we solve for the algebraic fields (one vector and three 
scalars) and reduce further over the gauge degrees of freedom (again one free function in 
the vector sector and three in the scalar). We are thus left with 3 types of gauge-invariant 
master fields, one type in each sector [36-42]; we call them h LaJ , f)^, 1 )(q. These appear 
in the action in similar forms (see appendix C, in particular (C.5, C.16, C.26)): 

4 (d+ 1 )(^- 1 )/ 4 (/- 1 )( <+ d) 

d(£+d+l)(£+d) (£+d+1)1 


N, 


Ld 


Seh = =- 

4^8(d+l)G d 


' dr r 


2t+d+\ 


+ d 2 c s (c s -d) f)io£f)Kri , 


(2.5) 


£ = 



21 + d + 1 
r 



( 2 . 6 ) 
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where N £ ( j = 9t Y^^+d /2 ) = ~( 2 t+d)\\ ' no ^ e that f° r d = 4 the scalar and vector expressions 
match correspondingly eq. (3.80) of [1]; there is no tensor sector in 4 d. 

For the inhomogenous (source) part of the action (2.3), we similarly decompose the 
sources T> 1V to T tt , T tr , T rr , T*, T r , T s , T s , T m , T rH , T K , as in (D.l). Collecting 
the sources into the combinations T, T matching the gauge-invariant master fields 
(see appendix D, in particular (D.15-D.18)), we find the action for the interaction of the 
fields with their sources is given by (D.14), and the complete action in our gauge-invariant 
variables is comprised of three parts 

S = S s + S V + S T , (2.7) 


where 


th 

for all sectors e 6 {5, V, T}, and 


S '=2 


N, 


G r iR: 


r 2i+d +^ _ ^* T e + c . c .) 




( 2 . 8 ) 


g d(£-\-d-\-l){£-\-d) y (d+1) {£-\-d-\-G)£ y 8(d+l) 

(Z^IK ’ 2(£-l)(*+d) ’ d^c s (c s -d) ‘ 


Variation of (2.8) yields the master wave equation (compare (2.14, 3.18, 3.26) of [3], and 
(2.9, 3.8, 3.49, 3.57) of [l] 6 ), 


0 = 


5S 


N, 


Si) 


Luj* 


G d R 


£,d r 2£+d+l 


2 , pfi , 2£ + d + 1 

uj + cL H-<9 r 


- n 


Luj 


( 2 . 10 ) 




2.2 Feynman Rules - Propagators and Vertices 

The solutions to the homogenous part of such wave equations are composed using the 
origin-normalized Bessel functions j^ +d / 2 (. ujr ) and ^ + y 2 ( ujr ) ( see a PP en dix E and [3]). 
Hence in the language of Feynman diagrams, we can represent the retarded propagator for 
the various waves using the diagram 



= G vet( r ', r ) = -iu 2 + G d CG j e+ j /2 (un) hf + ~ /2 (ur 2 ) ; (2.11) 


with n := minjr^r}, r 2 := max{r ; ,r}, G \ d = M f , ^ I?* ~ and 


M. 3 = 


IT 


e,d 2^+ 1 +^r(l + d/2) Y(£ + 1 + d/2) ' 


( 2 . 12 ) 


In particular j = 


dl\(2£ + d)\l 


i -i 


in odd d, M (J = § 


d\\(2£ + d)U 


i -i 


in even d. 


6 Notice that our subscripts i,d in R * j carry a different meaning than the subscript s in [l]’s Rl. The 
s in [1] denotes the field type, which for gravity is always 2; we choose to emphasize the 7, d dependence. 
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The vertices (source multipole terms) in the radiation zone are found by matching with 
the system zone according to the diagrammatic definition 



(2.13) 


where the blob on the right hand side means one should sum over all possible near-zone 
diagrams with an outgoing radiation leg. From the radiation zone’s point of view the 
sources Q LiJ are located at the origin r = 0, hence the radiation zone field can be written 
as 


f) 


e,EFT 

Llj 


(r) 



-Q 


Luj 


— VjJ 


2 i+d 


G d C\ 


£,d 


K+d^ ■ 


(2.14) 


On the other hand, in the full theory we may also use spherical waves to obtain the field 
outside the source as 


f)L0) 


/'dr'Tl UJ (r')G e ret (r',r ) 


(2.15) 

Jdr , 7L U (i J )j t+ j /2 (uj J ) 




From the comparison between (2.14,2.16) we read off the source multipoles as 

Ql M = J‘drj e+i/2 (u r )72 u (r). (2.16) 

Following [1, 26] we call this matching process a “zoom out balayage” ' of the original 
source distribution T^(r) out to the multipole carried out through propagation with 
Ji+d/v ( ujr ) • Thus we find both the source vertex and vertex for the doubled (hatted) source, 



(2.17) 


At this point we first introduce field doubling and employ the classically adapted Closed- 
Time-Path formalism; we work in the Keldysh basis. For a detailed description, see [1, 33]. 

It is worthy to note that the series expansions of both j i+ j/ 2 (E.3) and of the sources 
T[ u (D.15,D.17,D.18) include only integer powers of ui. This implies analyticity of the 

7 French for “sweeping away”; Poincare coined the term ([43, 44]), describing the process where a charge 
distribution in some spatial region is “swept away” to the boundary of that region, leaving the potential 
outside unchanged. 
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multipoles Q(uj) as functions of (complex) frequency, for any £ and d, which in turn results 
in local-in-time expressions for the time-domain multipoles Q{t ) - that is, the multipoles 
at a certain time t themselves depend on the energy momentum tensor at the same t alone 
and do not contain tails. The effective action S can still contain tails - this will happen 
in the case of non-even dimension [45-47]. Using the inverse transformations (D.2-D.5), 
source terms (D.15-D.18), integration by parts and the modified Bessel equation (E.2), we 
find the radiation source multipoles in the frequency domain 


Qs u 


Qkn 


{d+1)£ fd D xx L STF 


dc s (£+d+l ). 


H tt 



\d + 1 

+2iu;r [l 

+T a a 

f c s 

U+l 


-£+d- 


u 2 r 2 


+ rd r 


+£+d + rd r ) - oj 2 r 2 T rr 


2e 


(D) 

ttabke 


d D x 


£+d +1 

£(£-1) (D) ID) 


2 e 'Rabk l e ^a'b'k l i 


~h+d/2^r) , 

r . {e+d ) x bL-i T u dr (r e+d+1 j e+d/2 (u r )') +iu J T^x bcL - 1 ~j e+d/2 (ur) 

(2.18) 


1 d °x x WL 2 TZ a 'j e+dl2 (wr) 


Fourier transforming gives the multipoles in the time domain, recorded explicitly in (4.1, 
4.2, 4.3). These may also be written using the d-dimensional generating time-weighted 
function 8 f . (z) defined in [3] (following the 4 d definition [48, 49]), as 


Qe = 


d X X STp 


1-1 


dz8 { .( z)T e {f, u + zr) 


(2.19) 


The multipoles represent the gravitational source scalar, vector and tensor multipole mo¬ 
ments. A simple test of substituting d = 4 shows these multipoles coincide (for any £) with 
the scalar and vector multipoles given by [1, 26] (see (3.88) and (3.89) of [1]); there is no 
tensor sector in 4 d. In general d, for a system of N masses at leading PN order (low veloc¬ 
ities, weak potential), the leading multipoles (of order £ = 2, or quadrupole) of each sector 
are the mass (scalar, or electric) quadrupole, the gravitational source-current (vector, or 
magnetic, +1PN) quadrupole 8 , and the mass-tensor quadrupole (+2PN), respectively: 



i> m( 

A=1 N 

x l x^ — , 

(2.20) 


n 

f . . , \TF 


q%-> 

2 y ^m A 

(x l e J Hab x a v j A , 

(2.21) 


A= 1 




n 

- —, / 

aTF 


Qnn 

J2 mA { 

4 cd x c v) A ■ 

(2.22) 


A=i 


We note again at this point that in order to go beyond the +1PN approximation more 
Feynman rules, that correspond to radiation zone interactions, need to be included; see 
discussion in sections 2.1 and 3.3. 

8 The factor of 2 between the gravitational source current and the mass current, follows app. A.3 of [1], 
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2.3 Source Non-linearities 


In order to describe the radiative multipoles of a system of masses beyond the leading 
orders, we prescribe how the system’s T^ v is given by the masses and their coordinates. 
From the inhomogenous coupling (2.3), T^ v is the source for h /iv , represented as a 1-pt 
diagram (compare [1, 50]), 


T^{x) := -2 


Meff _ 

Shfj, u (x) 



(2.23) 


In the system zone we describe the metric using the three NRG fields -f 
defined by 



2A ■ dx 


_ e 2 <j>/d _|_ a ,^ d x l (lx^. 





(2.24) 


This dehnition almost matches that of [11, 35], but includes a factor of 2 in the dehntion 
of the coupling of the gravitomagentic vector potential A, as introduced in the 4 d case in 
eq. (A10) of [1]. The corresponding sources are given by T^ v i—> (p<p, J, £ j, defined by 


Smat — Prj> (f> J ' A 


-X ij a„ 


ij ■ 


(2.25) 


This implies 



(2.26) 


We refer to p$ as the gravitational mass density, to J as the gravitational source current, 
and to T, lJ the gravitational stress; together they comprise the NRG source fields. We 
further mark £ = £ aa , and note that T u = p^ + i£. In terms of these source fields we 
re-write the expression for the radiation source multipoles (4.1,4.2,4.3) as follows 


Qs 


(d+1) 


d(e+d)(£+d+ 1 ). 


' (I X X^j 1 p 


d (^±^±^±11 + rd}j~j^ /2 {ird t )p^ 
— 2r~ (yMdA ; 


x-J 


1 - 


+3e + d/ 2 ( ird t) ( r' ( /O 0 + ) - x a x b T, 


d 


\-\ab 


Qk 


e (D) . 

Habk e jjD x 




i+d+l 

(D) (D) 


/ j.bL—1 


y*(r\-d 

jZ) ^bb'L—2x~\aa' / 


J a + h+d/2^ ird t) T ‘ aCx 


r-\ac^bcL— 1 


e ha>b'k e , d XX S “ a (?, t)j l+ i / 2 Mi) 


(2.27) 
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We remark that upon substitution of d = 4, equations (2.27) reproduce equations (3.105) 
of [1] (using (E.2)). To use the NRG source fields defined in (2.26), we expand the diagrams 
perturbatively in PN orders. We now focus on the IV-body problem of masses {Tn A }^ =1 
whose trajectories are described by {®A(i)}^ = i, and influenced by their gravitational in¬ 
teractions. For such systems the gravitational mass density at leading and next to leading 
(+1PN) orders is found by Taylor expanding (2.24) and varying w.r.t cj) (using (2.25)) to 
be 





m A 5{x-x A (t)) 

A =1 



d + 2 
2d 


v A(t) - 

B^A 


G d m B 


| x A (t)-x B 


(2.28) 


When other (not only gravitational) interactions are involved in maintaing the bound 
system (see 1), they also contribute at +1PN. Their potential energy is to be included in 
the nonlinear correction (2.28), similarly to its last term (or in its stead, if the gravitational 
potential is much weaker then the PN parameter, (f>gravity «t) 2 ~ Mother )• 

Non-linear contributions to the additional NRG source fields J and 5W are to be 
derived in similar fashion. As we here calculate only up to +1PN, the leading orders of J 
and 5W suffice: 


J = 2 ^2 m A v A S (x — x A ) , (2.29) 

A= 1 
n 

Yi lJ = ^2 rn A v A v\ 5 {x — x A ) . (2.30) 

A=1 


These of course match (3.107,3.108) of [1] for d=\. 


3 Application results 


Having set up the apparatus of Feynman rules (by sector e and by multipole order £, in 
(2.2)), we now diagramatically construct the RR effective action [1-3]. We expand order 
by order in the PN expansion, i.e. in the small parameter v 2 ~ , and recalling that 

ur ~ w 1 / 2 counts as half an order. A general contribution to this effective action is given 
by a vertex-propagator-vertex diagram, 



0 )qLu* = _ iu} 2£+d GdC e^ Q Lu>QLu* + c c ( 3-1 ) 


Below we calculate the RR effective action to the leading and next-to-leading orders; we 
then discuss higher order contributions. 
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3.1 Leading Order 

The only contribution at leading order (LO) is that of the scalar quadrupole (’’Mass 
Quadrupole”, s = 0, £ = 2, also called E2), 


nLO _ 
b eff,u ~ 



—ioj 


d +4 


G * C U Q 


L 2 CJ / £ )L2UJ 

S ^S* 


+ c.c. 


4 G d . ^+3)(d+2) , Q L s 2UJ Q L sT 

2 d+4 r(3 + d/2)T(l + d/2) 


+ c.c. (3.2) 


Following the procedure of [1] ((2.57-2.61) & Sec. III.C.l), we can immediately find from 
this expression the RR force on the system’s constituent masses, the effective potential, 
and the energy dissipated as radiation. In the remainder of this section we will focus, 
for clarity, on time-domain results for even spacetime dimensions d (odd d), where (3.2) 
reduces to 


qLO ,even 
b eff 


(- ) ! ¥ d(d + 3K2+d) G d Q l s 2 d? +4 Q L s 2 
K J 2 d\\ (d + 4)!! t VS 


(3.3) 


The LO RR force is given by 9 


-rpi,LO,even 
b SF 


5S 


LO,even 

eff 


5x l (t) 


d- 1-1 

(—) 2 G d 


(-Y 


d(d+3)(2+d) 5Q‘J gd+^Qij 
4d!!(d + 4)H Sx^t) t s ’ 


(3.4) 


and it can be represented by a potential which generalizes the Burke-Thorne potential [51] 
to any even dimension, 


Vi v T en = {-)*& G d mx i x j d(d+3 j (2+rf) d? +4 Q)j 
BT ' 4dH (d + 4)!! 4 S 


The LO dissipated energy is given by 

dE LO,even d(d+3)(2 + d) „ 

- - = U) -X-X- (j T d 

du 4d!! (d + 4)!! 

which in the time domain gives the dissipated power as 


QW 


< p. 


,LO,even „_d(d+3)(2 + d) / 


rad 


> = 


4d!!(d + 4)!! 


G d (d t * Q'£(t) 


(3.5) 


(3.6) 


(3.7) 


using Qfj (t) to LO as given in (2.20). We note that for d = 4, (3.4) and (3.5) reproduce 
precisely the well-known results for the Burke-Thorne potential and RR force, and that 
for any d, (3.6) reproduces precisely the generalized Quadrupole formula of Cardoso et al. 
(1.1) [24] 10 ; (3.7) matches Einstein’s Quadrupole formula in d = 4 [25], and generalizes it 
to any even dimension. 

9 Notice a factor of 4 = | from switching from multi-index notation L 2 to usual indices ij; see Sec. 1.2. 

10 In order to revert from Gd to G, see Sec. 1.2. 
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Eq. (3.2) gives the leading order RR effective action (and allows the other quantities 
to be equally derived from it) in any dimension. Below we will discuss the higher order 
contributions; we compute the +1PN correction, and give an outline of higher order contri¬ 
butions. Although in this section we focus on even d, our frequency domain results hold for 
odd d 11 as well. The main difference in odd dimensions appears when Fourier transforming 
back to the time domain, where due to non-analytic features (branch cuts) of the effective 
action as a function of (complex) frequency, the transformation introduces non-local (in 
time) tail terms in the effective action (see Sec. II.C Sz IV of [3] for a thorough discussion 
and section 4 for time-domain results). 


3.2 Next to Leading order (+1PN) 

At the next to leading post-Newtonian order (NLO, which is +1PN), four effects must be 
considered; we shall label them E3, E25 1 , E2I5 1 and M2, as we now detail. The scalar 
quadrupole is supplemented by the scalar octupole (s = 0, 1=3, a.k.a E3), found exactly as 
in (3.2), but with £ = 3. In addition, the scalar quadrupole itself becomes more complicated, 
as it must include +1PN corrections (specialized to the IV-body problem), 

n 

Qs 2 ->• Qs 2 + ^Qs 2 = £ X A + 6l Qs 2 ’ (3-8) 

A= 1 

where 5 1 Qg 2 includes five possible corrections: the +1PN non-linear corrections to the grav¬ 
itational mass ( 5 1 Qg 2 NLl ), the contribution of the gravitational source current (5 1 Q^ 2 J), 
the contribution of a retardation effect from expanding the Bessel function j to 
subleading order (^Q^b), a term with double time derivatives of p<p (S 1 Qf j d 2 ), and a 
term with double time derivatives from the derivative of the Bessel function {5 l Q I £ ) 2 db). 
Altogether for any l, d we find the first correction to be 


8 1 Q l s = 5 1 Q l s NLI + 5 1 Q l s J + 5 1 Q L s b + 5 1 Q L s d 2 + 5 l Q L s db 


n 


=E 


mx 


d+2 2 
— -v z 
2d 


-E 

B^A 


G d m B 


\X — Xb\ 


! 2 (d+l) 
d(£+d) 


iuj(x-v) 


(3.9) 


uj 2 r 2 


(d 2 + {£+A)d+A) 


2d(£+d){2£+d+2) 



and specifically for £ = 2 (the correction relevant to +1PN order), This exactly coincides 
with (3.109, 3.110) of [1] for d = 4, as obtained in [52] and reproduced in [8] within the 
EFT approach. At +1PN, The correction 5 1 Q itself can appear in either the source 
vertex Q s (E25 1 ) or in the hatted vertex Q s (E2I5 1 ). 

Finally, the +1PN action includes also the leading order vector quadrupole (s = l, £=2, 
a.k.a M2), which includes different Q^ 2U Q^ 2UJ terms. We sum over the H indices to define 
the bi-vector multipoles 


Qm ~ 


1 


£+d+l.j 


d D x 


? A ( W5 {^h+d^r) ) J '+ ji+d/2^ r y- 




(3.10) 


One may also formally discuss the RR effective action in non-integer d, see [3]. 
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We accordingly define Qm , and replace C^~ by 


,-iM = nV Dt(d+ 1 , 1 ) 

^ D e (d+ 1,0) 


(3.11) 


which accounts for the summation over the different H combinations using (B.5), following 
[3]. At +1PN order, for the vector contributions we can set j^ +d j 2 i ujr ) = 1 and S =0; the 
other terms only enter at +2PN and onwards. 

With these definitions, the four contributions combine to produce the radiation-reaction 
effective action to NLO, given by 


qNLO 


—iu 


d +4 


—iuj 


d+6 


GdC 2 d ( Qs 2U1 Qs 2U * + S'Qs^Qs 2 "* + Qs ,u 6 1 Qs au *') 
G d C[ d QlMl. - G d C^ d * + c.c. 


(3.12) 


3.3 Beyond +1PN 

Our method allows identification and enumeration of higher PN order corrections as well. 
While we do not calculate them explicitly here, it is convenient to see how they too fall 
under the general prescription. Effects such as higher-order multipoles {£), higher-spin (s) 
sectors, retardation effects or non-linearities in the source terms all introduce integer PN 
order corrections, and will thus only enter at +2PN. Radiation-zone non-linear corrections 
and the effects of spin (intrinsic angular momentum) enter at different PN orders depending 
on the dimension, as we describe below. 


+2PN 

At +2PN, our formalism captures easily all the relevant near-zone contributions. The 
tensor field first enters at this order via the tensor quadrupole (s = 2, £ = 2, a.k.a T2). 
The vector sector now also includes the vector octupole (s = 1, £ = 3, a.k.a M3), as well 
as +1PN corrections to the vector quadrupole (M25 1 ) and (M25 1 ), which arise from the 
contributions of the gravitational stress (5 1 Q|(, i E), the contribution of two time derivatives 
of the gravitational current (8 l Qj^d 2 ), a retardation effect {5 l Q\ n b), and a derivative of 
the Bessel function (5 1 Q\ 2 db). The scalar sector now goes up to the scalar hexadecapole 
(s = 0, £ = 4. a.k.a E4), while corrections to the lower scalar multipoles include five analog 
corrections to the octupole (^ 1 <5f 3 9 2 NLl, 5 1 Qf s J, 6 1 Qf 3 d 2 , 5 1 Qf, i b, 6 1 Qf i db) and several 
new corrections to the quadrupole. For example, the contribution of T2 is given by 

= -i^ +4 G d Cl~Q^Q^*. (3.13) 

Radiation zone corrections 

We have thus far taken the propagator to be linear in the radiation zone (2.11). However, 
there are non-linear corrections to the propagator, which may be included through system¬ 
atic inclusion of non-linear terms in the action (2.2); see also discussion in section 2.1. The 
first of these contributions is represented by fig. 2, and is interpreted as scattering of the 
outgoing waves off the background curvature generated by the entire system’s total mass 
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M (labeling the vertex in fig. 2). Its value includes a factor of the gravitational potential 
GdM/X d , where A is the typical wavelength for radiation. As A ~ u; -1 , the value of such 
contributions is suppressed by at least (ur) d GdM/r d ~ v d GdM/r d relative to the leading 
order. In PN terms, this is equivalent to +(l + f)PN order. In d = 4, it implies a +1.5PN 
contribution; it is suppressed even further in higher dimensions. For further discussion and 
calculations for the 4 d case, see [1, 8, 10, 53, 54]. In the high d limit, linearized gravity 
suffices. 



Figure 2. Correction to the radiation reaction due to non-linear interaction in the radiation-zone 


Spin Effects 

Depending on the system’s parameters and nature of its constituents, spin effects may also 
need to be considered. This adds complexity and introduces new interesting effects which 
are worth exploring also in higher dimensions. The important scales determining the PN 
order - in addition to the typical orbital separation R and typical orbital velocity v - are 
the typical size of the spinning body ro, its moment of inertia / ~ mr q, and its typical 
angular frequency ui s \ together these give the spin 

S ~ I uj ~ m rg w s . (3-14) 

Spins can couple to the orbital angular momentum 

L ~ mvR , (3.15) 


or to each other. Couplings of the first type (spin-orbit, S-O) appear in the action (at 
leading order and up to dimensionless factors composed of the masses) as 

Sso = #so • Idt-S-S • L ~ [dt Grr ^v r ^’ s , (3.16) 

J r d+2 J R d R 

while couplings of the second type (spin-spin, S-S) appear (again at LO and up to the 
masses) as 

Sss = #ss • fdt^ 2 S • 5 ~ fdt , (3.17) 

where in the above equations we omit the indices of spin and angular momentum tensors 
since we are interested only in orders of magnitude. In terms of the post-Newtonian 
parameter v 2 , we recall that is the leading order, and that 


ro 

R 


> v 2 ' d 


(3.18) 
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where similarity occurs in the case of highly compact objects such as black holes (see 
[55]). We focus on two interesting cases: co-rotation (oj s ~ co ~ and maximal spin 
(oj s ~ ^ ~ ^). In the case of co-rotation, we find that S-0 effects enter with a suppression 
of v 2 (tj) 2 corresponding to PN order 1 + 2/d, while S-S effects enter with a suppression 
of v 2 (^=l) , which means PN order 1 + 4/d. For maximal spin, S-0 coupling incurs a 
suppression of v ^ (0.5 + 1/d PN) and S-S a suppression of (^) 2 (2/d PN). Thus for co¬ 
rotation, spin effects only enter beyond +1PN in any dimension. For maximal rotation, 
while in 4 d S-0 enters at +1.5PN and S-S at +2PN (compare [4, 12, 13, 56-64]), as the 
dimension grows both effects become more and more important, with both effects at +1PN 
for d = 5, and both entering before +1PN at d > 5. Also for d > 5, spin-spin interactions 
become more dominant than spin-orbit. 


4 Summary of results 


In general spacetime dimension d the radiative held decouples, at the linear level, into 
three sectors - scalar, vector, and tensor - with respect to the Ov j sphere. The radiation 
is generated by the corresponding multipole moments of the source. The linearized-GR 
multipoles are given in the time domain by 


Qs - 


(d+1) 


d(£+d)(£+d+1 ). 


’d°xx a L 


STF 


ri 


((£+d)(£+d+1 ) 


d + 1 


+ rd r + ^d 2 )jz + d / 2 (irdt)T tt 


—2 x a (l+d+1 + rd r 


’t)ut 


rpta 


Qxn 


+ 


((£+d)(£+d+1 ) 




d+1 


+ rd r je + d/ 2 ( ird t) Ta 


+x a xbj l+ d /2 (ird t )d; 


{2rj-iab 


2e 


C D) r 

Habk e l ( f> x 


£+d+l.j 


J++1~A 


/ rgbL —1 


rj-ita _ 


r &\-d 


h+d/2^ ird t) d t TaCx ' 


ac^bcL—1 


£(£-!) (D) (D) 

2 t Habk( t '2a'b / k i i 


ld D xx bb ' L 2 T aa 'j e+ d /2 (ird t ) , 


(4.1) 

(4.2) 

(4.3) 


where T^ v (r,t) is the energy momentum tensor. It is useful to sum over the H coordinates 
to replace the vector multipoles by 


Qm ~ 


1 


e+d+ij 


d D x 


f A ( ^ ( r ^e+d/2( ird t)) J'+ je+d^rdt)?- 




(4.4) 

These multipoles share the same form in any spacetime dimension. The Feynman prop¬ 
agator (2.11) introduces 2£ + d factors of ui into the RR effective action (3.1). In even 
dimensions, these are transformed to the time domain to become 2£ + d time derivatives. 

For the N-body problem at leading post-Newtonian order, the radiation-reation effec¬ 
tive action includes only the leading order of the scalar quadrupole Q 2 g = Xm=i m A(x 2 ^) TF 
(2.20), and the action reproduces the generalized Burke-Thorne potential, as well as eq. 
(1.1), as described in Sec. 3.1. 
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At +1PN order, contributions arise from the scalar (mass) octupole, 


Qs 3 = Qs k = y ~]mA 


A= 1 


x i x j x k - 


1 


D + 2 


'S - ^s 

from the vector (current) quadrupole, 

n L 2 — O— 2 

Vm — Vm — z 
and from corrections to the scalar quadrupole, 

Gd m B 


6 ij x k + + S jk x l )x 2 


J A 


m [ m ( FA ^ (v) 

A=1 


J A 


(4.5) 


(4.6) 


A=1 


d+2 2 
—- Va 
2d 


£ 

B^A 


\xa-xb\ 


,x *- 2 $£) a, ( 3A -* AX *) 


+ 


d 2 +6d+4 


2d(d+2)(d+6) 


d't (taXa 2 ) 


(4.7) 


Taken together we find the radiation-reaction effective action to +1PN, which in the case 
of even dimension d is 


S? f ? = (-)^G d dt 


( Qs 2 d? +1 Q L S 2 + Qs i dt +1 6 1 Qs 2 + 6 ] Q^d? +l Q 


d{d + 2)(d + 3) 

2d!!(d + 4)!! 

d(d + 4)(d + 3) ~l 3n d+3 n L 3 . 2d(d + 3) f ) L 2F) d+i n L 2 
6 dl!(rf + 6)M Qs ‘ Qs + 3dU(i + 4)U QM '‘ Qm 


(4.8) 


This is the main result of this paper. It matches (3.111) of [1] in d = 4 (upon re-introduction 
of l/£\ from the summation convention), and extends it to different dimensions. The energy 
dissipated from the system, as well as the gravitational self-force acting on it, can be read 
off from this equation. 

While all the results of this paper, presented in the frequency domain, are valid for any 
spacetime dimension, for odd dimensions the Fourier transformation to the time domain 
introduces non-local “tail” expressions in the action. This happens due to the appearance 
of branch cuts in the frequency domain effective action. For a thorough discussion of this, 
see eq. (4.7, 4.8, 4.14, 4.15) of [3]; as an example, the time-domain effective action in odd 
d is given by 


Seff = G d dtj^i 


\<+ 


d +1 


C s l4 S& (t) + C% SM (i) + S^ (t) 


(4.9) 


S e (t ) = Qi(t) 


\H(2l + d)-H{l+±)\d? + «Q L e {t) 


- f dt 1 

J —OO 


1 


t - f 


d 2 J +d Q k (t’) 


(4.10) 


regularized 

where e G {S,V,T}, H(n) is the n’th Harmonic Number, and regularization (hence the 
effective coefficient of the local part) depends on shot-distance details of the system, in 
particular on the scale (see also [3, 19, 65, 66]). 
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A Properties of spherical derivatives 

The generalized spherical metric for coordinates {t, r, Oi, P 2 , ■ • •, P^ +1 } is 


i— 1 


g,,u = diag{-1 , l,r 2 n 1 ,r 2 n 2 ,... ,r 2 n^ +1 }, U,; = JJsin 2 ^-, 

3 = 1 

from which we find all the Christoffel symbols 


n a n. 


^ S^a _ ^ rfta T't _ "pCt 


r Q^l a ^b’> rQh ^ rQ a ™ ^ olB t/3 


/r> u u O /u 


= 5^ a 9(a — c)cotQ c + 6^ a 9(a — b)cotflf> — 6^ b 9(b — a)—- cotH a , 


U 

n„ 


where 9(a — b) = 1 for a >b, otherwise 0. The Riemann tensor on the sphere is 


(A.l) 


(A-2) 


dio. a fl b ^l c fl d 9^ia^c9Cl b fl d 9£l a Qd9Qb^c • 


(A.3) 


We shall also make use of the relation 

g nn 9 W = d + 1, (A.4) 

and note that 

gw Si = 0 . (A.5) 

The volume element is 

d 

\J—g d d x = r bl+1 (sin Pj) d+1_ * dtdrdVtv • •f2^ +1 = dtr^ +l dr dQ^ +1 . (A.6) 

2—1 

Using the Riemann tensor we find the commutation relations for a tensor on the sphere: 

[Dn a ,D n J Vq c = SUhsSi = 9n c [n a Vn b ] (A.7) 

The commutators for tensors give 

\Du a i 9Q a Q c TQ b ci d — gQ b ci c Tci a n d + gn a ci d Tfi b Q c gci b ci d Tci a ci c , (A.8) 

[Dq". Dq\ (DqhVp) = ((1 + 1)DqVp - gQpDn»Vn'i . (A.9) 
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B Spherical fields 


A-olP 


=£ 


A t dan L +A m n^ n 


For a rank-2 tensor A^ u in d dimensions we use the following spherical decomposition (A’s 
L, lo indices are suppressed): 

^ Att A tr A t n 

Ai r A 

rr A r n 

\A-tn A r n Aqq/ 

(A tt n L A tr n L 

ArrTlL A r del'll p, + A r ^fl^^ 

As'> l fi£i'+AsnnQ,+Axnx£ l w+An'2n^ nn , 

where we use the scalar multipoles n L , the divergenceless vector multipoles n ([+ 12 and the 
tensor multipoles n^/nn' (which are symmetric, traceless and divergenceless) 13 . These 
multipoles are all dimensionless, and depend only on the angular coordinates. They are 
related to the scalar, vector and tensor spherical harmonics: 

, r ,L 


—iut 


(B.l) 


V 


n = —r = Y , D n n = Y n , n nn , = g nn >n , n an , = D n D n * + 


d +1 


9w n , 


L _ v 

n m — ~r ~ 


i > n KQf 2 ' — 9 + Ar^n) > n H3nn'~ ~ ■ ( B -4) 


Using [3, 72, 73], the numbers Dt(d+ l,e) of such independent multipoles of order £ and 
sector e on the O Vj-sphere are 


De(d + 1,0) — 


{2£ + d)(£ + d - 1)\ 


Di(d + 1 , 1 ) — 


£{£ + d)(2£ + d)(£ + d - 2 )! 


d\£\ ' " ' ' ' (d- 1 )!(£+ 1 )! 

(d + 2){d-l)(£ + d+l)(£-l)(2£ + d)(£ + d-2)\ 


Di(d + 1 , 2 ) — 

2d\(£+ 1 )! 

The spherical multipoles satisfy the derivative properties (using appendix A) 
dthLu — iuih , drhLu — ^Lui > d r Tiyr — ® i drXx — ~Xx 


DfihLu — 0 , — 0 , Dfin^nn 1 — — 0 


CIQ 1 L _ « 

9 n nnn r ~ u 


j Dnnhw — Dq> n L , Dnn^ nn , — ——n#Q, 

d+1 z 


(B.5) 

(B. 6 ) 

(B.7) 

(B. 8 ) 


2 The vector multipoles are enumerated by an antisymmetric multi-index H taken from the Hebrew 


alphabet, representing D — 3 spherical indices: 


L (d+1) r\Cl r L 

- ~ , D n = 




(*(r 


AV) 


n 
/h ci 


(B.2) 


where . is the completely antisymmetric symbol on the fij +1 -sphere, A is the exterior product and 

* is the Hodge duality operator [67-70]. The spatial Levi-Civita tensor will be marked ei+.„ n ■ 

13 The tensor multipoles (of rank 2, generalizing [71]) are enumerated by 2 antisymmetric multi-indices 
N,N': 

L (d+1) (d+1) t~\ + t~v \F ' L /t~) Q \ 

2f2' 'Cl'ty'D D 'Cl • (B.3) 
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where a prime denotes r-derivatives (' ■= d r ). They are eigenfunctions of the Laplace- 
Beltrami operator on the (d + l)-sphere, A ^ +1 = DqD^, with eigenvalues 

A d ' + 1 n L = -c s n L , A d - + 1 n^ n , = -c s n^ n , , A j +1 D n x L = -( c s -d)D n n L (B.9) 
A d+i n ^n = -(c a -l)n£ n , Aj + 1 n£ nn , =-(c a -d)ra£ nn , , (B.10) 

A d+i™nn' = -(c s -2(d+l))nm' > = _ ( c s _2 ) n K3oo'- ( B -ll) 


The scalar, vector and tensor basis elements are orthogonal to each other, and we use the 
following normalization conditions in d dimensions ([72, 73]), 


Jd^d + i nL e(^d+i) n l ' (^d+i) = ’ 

jd^d+\3 Do,nL e Dfi'n >’ = c s ■ N f ,jQ ( j +1 8fji5 LtL ' fl > 

Jd£lj +1 g n nh n wn' = Cs ' ^td^d+i > 

K, 5 = (d+ l)A^jllj +1 S((/S LeL i fi , 


.nn'_ dc s c s 


'd^ +l9 ^g^n- 

SlQ.''Sr9'L t L 


^e^d+idii'^LtL',, j 




y _ C s c s 

n'S> n K , n i 'S>' ~ ~A~ 


^ | ^ t,a a-(-i ^ 

' ^e,d^d+l^ U 'dL e L', j 


dn ^fi^ nn, ^’* , ’ Sr,n ran« n Na'n'*' = d 2 c s (c s -d)-N e fL d+l 8w5 Ll L' <W<to' ■ (B.12) 




B.l More spherical expressions 


The d-dimensional trace of ^Lj, is given by 


A:=g aa A a „=y' I -A„ + A rr + ) n L e"^ ^TAr^e**, (B.13) 




and its divergence is given by 

(dhM) a = A a pf =—dtAt a +{A ra j—Y^ a A lxr )-\— _ r^/^4 /iQ; ).(B.14) 

Explicitly, its components are given by 

(di vA) t =^e-^ 


■ a I n d + 1 \ C s . 

tujAff + \ d r -\ - ]Atr - kAi 

1 r / r z 


til, 


(B.15) 


(diw4) r = 


—iujt 


( _ cZ H- 1 \ Cg . d - (-I. 

iujAt r + I d r H-I A rr - ^A r -o— A§ 

1 r r z r 3 


n L , (B.16) 


(divA) n = 


—iu)t , 


w^+(a r +— dCs A S 

l r r 2 (d+ 1 ) r 2 


Dlitil 


+ 


. / d +1 \ c s . 

iujA t K + I o r H—— I yr r K — ^2 ^ 




(B.17) 


19 




















where we have used 


D Aw = 


A s - Dqu 1 " - -^n^iAx 


(B.18) 


C Action in spherical fields: Kinetic terms 

In this appendix we express the linearized GR kinetic term with gauge-invariant spherical 
fields. We start with the action (2.2). Note that the terms containing traces of the 
perturbation ~ Dh% can be expressed only with d-scalars, and thus do not contribute 
to any vectorial or tensorial terms; and that in fact tensorial terms arise only from the 
term ~ D^h a ^ Dryh a a. After expansions, partial integrations and use of the normalization 
conditions, we find the homogenous part of the action (L, u indices will be suppressed 
throughout much of this appendix), 


Seh = 


JdrL(r) = 


N, 


%Z{d+l)G d 


d+1 


drL(r ) 


L]uj L,uj 

where the action decomposes into scalar, vector and tensor sectors 

L(r) = L s (r ) + L v (r) + L T (r). 

The tensor sector includes a single contribution, 


Lt = 

d 2 c s (c s -d) 


Enn = 


2 r 4 


u/+af+ 


2 , d —3 c s + 2 (d— 2 ) 

” (Jr 


r 

~-2 1 


hxn- 


This matches (5.6) of [27], if we set ~ r 2 h *q. Defining 

= r~ ie+2) hn2, 

it can also be presented in the master-form 


Lt 


X 




8(d+l)Gd 


where £ is the master operator defined in ( 2 . 6 ). 
The vector sector is comprised of three terms, 

Ly = h^Em + h*^E r n + h^En, 


Em = —9 


|d 2 + — d r - Cs+d n \\dr + :1 - 1 1 


d -1 




d+1 \ 

r J 


2 r 2 


En — 


r 

CsC 




Cs 

2r 2 


S'-'S 


4 r 4 


2 (d r +— 1 

h r ^ 

2 d—3 . „ 2 —d 

lo -\~d r H- d r +2 — 

h'K 

r \ r J 


rjrt J ‘ X. 



(C.l) 


(C.2) 


(C.3) 


(C.4) 


(C.5) 


(C.6) 

(C.7) 

(C.8) 

(C.9) 
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The field h r $ appears in the equation E r $ = 0 (C. 8 ) without derivatives. Thus we can 
solve for it algebraically, 


hrH — 


My(r) 


Vjjr 


d r --)h m - <9 r--)hn 


(C.10) 


where My(r) = a; 2 r 2 — c s (compare (5.9, 5.10) of [27]). We thus see h r # is an auxiliary 
field. We plug this expression in the equations for E t n, E#, Ly, 


L v = h* m E m +h* H E H 


Co Co 


4 r 2 My 


(r^)*2y (r 2 <f> K ) = 


Co Co 




r d+3 




,d +1 


-d r 


My 


dr®H 


(C.ll) 


where = <9 2 +(d—3 — jffyJ l<9 r +(2— + and the vector part of the action is 

seen to depend only on the single gauge-invariant field (compare (5.12, 5.13) of [27], and 

[41]) 

= r 2 (2 htn+iuhn). (C.12) 

Since Ly only appears in the integral (C.l), we can integrate by parts to obtain 

„2 


Ly = 


Co Co 


($*) ~ 


My(r) 


(d r ^y 


Defining the field using the canonical transformations 

dL c,c 


-l)(e+d)r i+d+1 \) H y.~- 


-1 )(£+d) / +d+ 1 tiK) , : = 


■'S^s d+ 3 




a I 


2Ady 


(*«)' 


4 tt 






= 1 c s c s r d+1 $ 
2 

=> = 


e(i+d+1 ) V r 

K 

4 / t+d+l 


£(£+d- 1 - 1 ) r d+1 


fm) , 


(C.13) 


(C.14) 


(C.15) 


and again using integration by parts, we can present the vector contribution as well in the 
master-form (compare with (C.5), differing only in the prefactor constants) 


= 8 ( 1 — !) (£+d) r 2e +< i+i ^ 

8(d+l)G d (£+d+1)£ 


(C.16) 


In the remaining scalar sector, we also expect a reduction to a single gauge-invariant 
field. We first notice that 3 of the fields (hf r , h rr , h r ) appear without r-derivatives in 
the action; thus variations with respect to their complex conjugates (hl r , h* r , h*) give us 
algebraic equations 11 . Hence these fields can be found by solving the system given by the 

14 The complete derivation of the scalar sector of the action, along with the collection of terms to the 
gauge invariant combination, can be found in Mathematica code in our additional notes online [74]. 
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(hermitian) matrix A 15 , 



^ hfr \ 

( Cs 

— (d + l)iuir iujr \ 

(Bi\ 

A 

hrr 1 — ^4 .— 

(d + 1 )iur 

d{d + l) d , B - 

b 2 


Vfhr) 

\ —icor 

d uAr 2 +2d J 

Cs ' 

W 


with 


B\ := c s d r ht — (d+l)zw —jhs 

c s -(d+l)rd r . 

02 ■■= --- htt+iuc s ht - 


fd(c s +d+l) 

d(d+ 1) 

1 2 r 2 

2r 


B 3 := (rd r - l)h t t+ico(rd r - 2)h t + -(- - d r )hs-rd r 

r \ r 


We solve these simultaneously to find 

^ htr 


d r -— u 2 Jhs — 

M • 


dc. 


(d+l)r 2 


dc s c s hs 

2(d+l)r 2 

(C.18) 


h rr I = A 1 B = 


A~ 2 B, 


2 M s (r) 


(C.19) 


i - 1 = 


^d (2dc s — (d+l)u; 2 r 2 ^ 2iur(d(c s — 2d— 2) — (d+\)uj 2 r‘ 
* —Ac s d 


V 


2 c s 


— iurc s d{d+l) 


uj 2 r 2 — dc. 


f= * c s (d+l)^2(d+l)ij0 2 r 2 — c q d 

Ms(r ) := c s c s d 2 — 2c s d(d+l)oj 2 r 2 + (d+l) 2 u; 4 r 4 , 


\ 


where the matrix terms under the main diagonal (in asterisks) are completed by hermiticity. 
Plugging (C.19) back into the action (C.1-C.2), we find that hs drops out entirely, and the 
three remaining fields htt, ht and hs appear only combined in the single gauge invariant 
field (compare (4.30) of [27] in the flat space limit) 

$ := htt + 2 iuht — (j0 2 hs ■ (C.20) 


The scalar action is thus given by 

L s = —Vt + (3 (<9 r M+$*d r $) + 7 4>*$], (C.21) 

M s {r) 

a = — d(d+l)c s c s , j3 = —— (c s d— (d+l)w 2 r 2 ^ , 7 = —^ ^2(d+l)w 2 ?’ 2 — c s dtj . (C.22) 


Following [27] we canonically transform to 4? with the generating function, 


F 




-c,r 2 




(C.23) 


15 We mark that this matrix is identical to (4.35) of [27], adapted to a flat background - up to ordering 
of the fields, scalings in the field definitions, and sign mismatches stemming from the choice of Lorentzian 
over Euclidean signature. Using our definitions, the matrices A, A -1 , A _1 are dimensionless. 
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and after partial integrations find 


L S = _ d(d+l)c s ^ dt+u , + [cs + 


8(d-\-l)Gd c s 


d 2 - l\ 1 


4>, 


(C.24) 


which we recognize as the Zerilli action [38]. Defining the field 

(C.25) 

we find the homogenous scalar action in the familiar master form (compare (C.5, C.16)) 

N t,d 4 (d+1) (£-!)£ 


= d(£+d+ 1) 

; 2£ 


r 2W +V^ . 


(C.26) 


Ls = 

8(d-\-l)Gd d(£+d+l)(£+d) 

D Action in spherical fields: Source terms 

We wish to construct the inhomogeneous part of the action S ma t (2.3), which describes the 
interaction of the gravitational field with matter sources, in terms of the gauge invariant 
spherical fields and corresponding source functions. We decompose the stress-energy tensor 
as 

T tt =^Ttn L e-^ , T tr = ^ 

T m = ^ e~ iu}t , T rU = ^ ( T r Lu] d n n L + T^n^) e~ iujt 






TlJf + ftn™ 1 ' + Ttn« L nn ' + 


—iut 


(D.l) 


We shall also use the inverse transformations, 


-,(tt/tr/rr) , s _ 


' Luj 


d^d+i 


dn- M dte iut n L T Wtr/rr \r,t), 


M = K 


d +1 


T% r) *(r) = 

^(r) = 


(D.2) 

(D.3) 

(D.4) 


1 


idte^ei^Le^L, l -^J^n bb ' L - 2 T aa \ (D.5) 


N f &;,J 2 c s {c s -d)JJ d+1 * abke ^ a ' Vktf r 2 


£,d d+1 

where T^ r = T r/i = and T rr = ^T cd . 
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Using the divergence expressions (B.14-B.17), current conservation D^T^ V = 0 is recast 
as the d equations (we henceforth suppress L, ui indices) 


0 = -iuT tt + ( d r + d H± jT tr - c s T t , 
d + 1 


0 = -iuT tr + | d r + 


, r -r —- I T rr - c s T r - (d + 1 )rU 


0 = -iuT* + f d r + — ) T r + T 3 - 4^f s , 

V r J (d+ 1) 


0 = -iuT m + ( d r + — ] T rK - -T k 


(D.6) 

(D.7) 

(D.8) 

(D.9) 


From these we can eliminate the fields T t ,T r ,T s , and T**, retaining only the scalars 
T“,T tr ,T rr ,T 5 , the vectors T m ,T rK and the tensor T ra : 


rj-it _ _ 

C s 

1 


-iuT u + I d r + |T* r 


T r = 


rpS _ 


(d+1) 

d Co Co 


-iu>T tr + [ d r + d — ] T rr - (d + 1 )r T s 


(D.10) 

(D.ll) 
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i tr 
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+ ^ + 2 — ^ ' T rr 


rjpN _ 


-iuT m + I d r + 


d+3\ 

r I 


-irH 


(D.12) 

(D.13) 


Plugging the expansions (2.4,D.l) into (2.3) and using (D.6-D.9) and (C.10,C.19), we 


find 


Sm.nt — - 


HI 


t)T*+ t)*T**+ t)iaT^*+ c.c. 


*a 


dr , 


(D.14) 
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where the source functions (Scalar, Vector, Tensor) are given by 

N e ,^d + i (d + l )/+ d + 3 


T := — 


d(£+d)(£+d+ 1) 


d o d[ c s 2{d+l 
—d r -\- uj — 


(d + l)r 2 
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(D.15) 


J := T tt + 
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(D.16) 

(D.17) 

(D.18) 


E Origin-normalized Bessel functions 

For B a Bessel function of the first or second kind or a Hankel function, i.e. B 6 {J, Y, H ^j, 
and with a representing its order, we define the origin normalized Bessel functions b a as 

t B a (x) 


b a := T(a+ l)2 c 


x v 


These functions satisfy the equation (compare 2.10) 

b a (x) = 0 . 


•-)2 i 2a + 1 

dj, H- d x + 1 


(E.l) 


(E.2) 


The purpose of the definition is to have j a normalized to 1 at the origin x = 0. Around 
the origin it is given around by the Taylor expansion (which contains only even powers) 


x 


^(2p)!!(2p + 2a)!!‘ 


+ ••• • 


(E.3) 


2(2a + 2) 

The asymptotic form for x —> oo is best stated in terms of the Hankel functions h± := j ± iy 

(E.4) 

Y M X~ 1 - 

For more details, see appendix B.2 of [3]. 
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